Many-body effects on the Pxx ringlike structures in two-subband wells 
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The longitudinal resistivity pxx of two-dimensional electron gases formed in wells with two sub- 
bands displays ringlike structures when plotted in a density-magnetic- field diagram, due to the 
crossings of spin-split Landau levels (LLs) from distinct subbands. Using spin density functional 
theory and linear response, we investigate the shape and spin polarization of these structures as a 
function of temperature and magnetic-field tilt angle. We find that (i) some of the rings "break" at 
sufficiently low temperatures due to a quantum Hall ferromagnetic phase transition, thus exhibiting 
a high degree of spin polarization 50%) within, consistent with the NMR data of Zhang et al. 
[Phys. Rev. Lett. 98, 246802 (2007)], and (ii) for increasing tilting angles the interplay between the 
anticrossings due to inter-LL couplings and the exchange-correlation (XC) effects leads to a collapse 
of the rings at some critical angle ^c, in agreement with the data of Guo et al. [Phys. Rev. B 78, 
233305 (2008)]. 
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The fascinating quantum Hall regime hosts a number 
of fundamental physical phenomena, being also relevant 
for metrology (standard for resistance) and as an alter- 
nate means to precisely determine the fine structure con- 
stant. The spectrum of two-dimensional electron gases 
(2DEGs) in the quantum Hall regime is quantized into 
highly degenerate Landau levels (LLs) pLj. At opposite- 
spin LL crossings near the Fermi level, a ferromagnetic 
instability of the 2DEG may arise thus leading to a quan- 
tum Hall ferromagnetic phase. This spontaneous spin po- 
larization of the electrons lowers the repulsive Coulomb 
energy of the Fermi sea, because electrons with parallel 
spins avoid each other due the Pauli exclusion princi- 
ple. Even for vanishingly small Zeeman splittings the 
exchange-energy gain can stabilize quantum Hall ferro- 
magnetism at low enough temperatures [2]-[4]. 

Quantum Hall ferromagnetism has been extensively 
studied in the quantum Hall regime via magnetotrans- 
port measurements [5^. Near LL crossings in tilted mag- 
netic fields 5, the longitudinal resistivity pxx vs B of 
wells with a singly occupied subband exhibits ubiquitous 
hysteretic spikes, which signals quantum Hall ferromag- 
netism [GHS]. In two-subband wells, spin-split LLs from 
distinct subbands cross even without a tilted B field and 
can form closed loops [ABCD loop, Fig.[TJa)]. Quite gen- 
erally, pxx is directly related to the energy spectrum near 
the Fermi level (linear response) and hence the topology 
in Fig. 1(a) translates into ringlike structures [H [10] in 
pxx when plotted in a density-5-field diagram 77-20 — 
Fig. [ijb). Recently Zhang et al. [TT] have shown that 
near opposite-spin LL crossings the rings "break" at low 
enough temperatures (70 mK). NMR measurements [12] 
near the broken edge C show a high degree of spin polar- 
ization, which points to a ferromagnetic instability of the 
2DEG. For tilted B fields some of rings "shrink", fully 
collapsing for angles above a critical value [T3] . 

Here we use spin density functional theory (SDFT) [T4] 



together with a linear response model [T5^ to investigate 
the shape and spin polarization of the ringlike structures 
in realistic quantum wells with two subbands at various 
filling factors i^, as a function of the temperature and tilt 
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FIG. 1: (a) Two-subband GaAs well and schematic fan dia- 
gram with LL crossings from distinct subbands. The ABCD 
loop gives rise to ringlike structures in the calculated niv,-B 
diagram of pxx (b) for v — \. This ring "breaks" at lower 
temperatures (c) due to quantum Hall ferromagnetic transi- 
tions, thus displaying a high spin polarization within (d) ^12] , 
For increasing B field tilt angles ^, the ring shrinks (i.e., A 
and C move closer) and fully collapses (A=C) at Qc (e), in 
agreement with the data ^ISj (cf. empty and solid circles). 
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angle of the B field. We find that exchange-correlation 
(XC) effects are crucial to quantitatively describe the ex- 
periments: in particular (i) the v = A ring breakup at 
low temperatures [Fig. [ijc)] follows from quantum Hall 
ferromagnetic phase transitions [10j[TT], Fig. 2. The cal- 
culated spin polarization [Fig. [ijd)] within the broken 
ring reaches 50%, consistent with NMR data [12 . (ii) 
The shrinkage of the v = A ring for increasing 6 and 
its full collapse = Oc [Fig. [ije)] [13] arise from the 
interplay between the anticrossings due to the inter-LL 
couplings and the exchange field, Fig. 3. We note that 
only rings formed from consecutive LLs, for which inter- 
LL coupling is operative, collapse for increasing 0. 

The quantum phase transitions we find here are not 
specific to the v = A ring. They are general and should 
also occur for = 6 and others, but for distinct ranges of 
parameters. Other 2DEG systems, e.g., formed in Mn- 
based wells [16], can also show peculiar ring structures. 

System. We consider the structures of Zhang et al. 
[TOHlS]: a wide 240 A GaAs square quantum well with 
Alo.3Gao.7As barriers and symmetric (5-doping (Si) with 
240 A spacers [Fig.jIJa)]. The electron density in the well 
is controlled by a gate voltage, as in an ideal capacitor 
model [17]. At zero bias n2D = 8-1 x 10^^ cm~^. The elec- 
tron mobility /ig = 4.1 x 10^ cm^/ Vs is assumed constant 
for the entire B field and gate voltage ranges. 

Kohn-Sham problem. The Kohn-Sham implemen- 
tation of density-functional theory maps the problem 
of fully interacting electrons onto a non-interacting 
Schroedinger equation - the KS equation - with electrons 
in an effective single-particle potential [14 . For magnetic 
fields B tilted with the 2DEG normal, this reads 
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where m (0.067mo) is the effective mass, ge (—0.44) 
the bulk g-factor, Px,y,z the x, z components of the 
electron momentum operator, ujc = eB cos 0/m the cy- 
clotron frequency, = ± (or t^i)^ = ePsin^/m, 
-i^Py/h^ — h/eBcos9 the magnetic length and 



v^ffiz) = Vc{z)^vh{z; [n]) +<^(z; [n^,nj), (5) 



Vc{z) is the structural well potential. The Hartree poten- 
tial vh{z; [n]) is obtained self-consistently from Poisson's 
equation. For the XC potential v^^{z; [n>^,n^]), we use 
the PW92 parametrization [18] of the local-spin-density 
approximation (LSDA) [19]. Here we have approximated 
the electron density n{x^y, z) by its average over the xy 



plane n{z) = n^{z) + ni{z) [8]. This renders both the 
Hartree and the XC potentials dependent upon only z. 

Perpendicular B field. For ^ = 0°, cjp = =^ 
SHq = 0, the KS equation ([T]) is separable in the xy 
and z variables and has eigenfunctions ip^^^ ^^{x^y^z) = 

eyiY>(\kyy)ipn{x)x1" {z) (Landau gauge), with (pn{x) 



being the n-th harmonic-oscillator eigenfunction centered 
at xq = —hky/muc and ky the electron wave number 
along the y axis; Ly is a normalizing length. The KS 
eigenenergies are = Sn -^^^^ + del^BC^zB /2 ("Landau 
fan diagram"), with = (n + 1/2) hwc^ n = 0, 1, ... the 
LL energies (degeneracy ub = eB/h) and e^^ the quan- 
tized levels obeying H^^Xi^ = ^i^Xi% ^ = 0, 1, .... with 
a self-consistently calculated chemical potential fi. 

Tilted B field. For 6> > 0° the KS equation ^ is not 
separable because SHq ~ sinOzPx ^ 0. However, since 
SHq <C i^o = ^11 only couples consecutive LLs 

from distinct subbands, we can obtain the KS solutions 
tlj{x^y^ z;0) as an expansion in terms of the eigenfunc- 
tions ^ (x, z; 9) of Hq. We perform this expansion 
at every iteration of our self-consistent scheme. We ob- 
tain good results by truncating the expansion for energies 
greater than /j. -\- ksT. This LL coupling leads to anti- 
crossings of the KS energies for equal-spin LLs, which 
ultimately make the ring shrink for tilted fields. Fig. 3. 

Linear-response pxx- By assuming that the KS 
eigenvalues represent the eigenenergies of the actual 
(Fermi-liquid) quasi-particles in our 2DEG, we use them 
in a Kubo-type formula [15 to calculate the conductivity 
tensor a. For instance, within the self-consistent Born- 
approximation with short-range scatterers [15 ~ 
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Text denotes the width of the extended-state region 
within the broadened density of states and f{e) the 
Fermi function. We obtain the resistivity from p = a~^. 

Spin-polarized rings. Figures [TJb) and [TJc) show our 
calculated n2D — B diagram of pxx for two different tem- 
peratures T = 340 mK and T = 70 mK, respectively, near 
the z/ = 4 ring. Similarly to the experiment of Ref. [11] , 
we find that the u = A ring "breaks" at the opposite spin 
LL crossings (points A and C) at lower temperatures [20] , 
Fig. [TJc). Figure [TJd) shows the corresponding n2D — B 



diagram of the spin-polarization ^ = (n, 



^2D = ^2D+^2D- Interestingly, we find a high spin polar- 
ization 50 %) within the u = A ring. This high polar- 
ization and, more importantly, its abrupt variation at the 
opposite-spin crossings A and C signal quantum Hall fer- 
romagnetic phase transitions. For T = 340 mK the spin 
polarization of the u = A ring (not shown), though high, 
varies smoothly at the opposite spin crossings. We note 
that the high spin polarization ^ within the ring points 
to quantum Hall ferromagnetism, being also consistent 
with resistively-detected NMR data available [12]; how- 
ever, we contend that the high ^ and the discontinuities 
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FIG. 2: (a) Longitudinal resistivity pxx calculated for n2D = 
7.3 X 10~^^ cm~^ [see horizontal dashed line in Figs. 1(b)- 
1(c)] for T = 70mK (dotted line) and T = 340 mK (solid 
line, slightly shifted upwards for clarity). Note the distinctive 
spike in the higher temperature pxx B ^ 7.6 T. (b) and (c) 
show the corresponding Landau fan diagrams for low and high 
temperatures, respectively. Note that for T = 70 mK, the LLs 
show discontinuities near B — 7.6 T. These discontinuities 
suppress the pxx spike at B ~ 7.6 T because the corresponding 
levels do not actually cross the chemical potential p - they 
suddenly jump over it. At higher temperatures the Landau 
diagram is continuous and the spike appears m pxx • 



of Pxx at the crossings A and C [Fig. 1(c)] constitute the 
signature for the quantum Hall ferromagnetic instability. 

The contrast between the low and high temperature 
results is more clearly seen in Fig. [2j which shows Pxx 
for n2D = 7.3 x 10^^ cm-^ at 340 mK and 70 mK. The 
spike near B = 7.6 T [see arrow in Fig. [2|a)] comes from 
the left edge of the u = 4 ring (point A in Fig. [T]) and 
is suppressed at T = 70 mK. The Landau fan diagrams 
for both temperatures differ substantially only around 
this region [see arrows in Figs. [2|b)-(c)]. At T = 70 mK 
the diagram shows an abrupt transition and the chemical 
potential p jumps to the spin-down state of the lower 
subband, thus suppressing the pxx spike. Note also the 
exchange enhancement of the spin splittings in Fig. 2(b)- 
(c) when p lies essentially between the spin-split LLs. 

A relevant parameter in our simulations is the LL 
broadening F. For short-range scatterers, the electron 
mobility pe and the LL broadening are related by F = 
TQ^B/pe [21 , with Fo = (2/7r)^/2^e/m. We use F70 = 
0.130V^ meV and F340 = 0.150V^ meV to simulate the 
ring structures at T = 70 mK and 340 mK, respectively 
[see Figs. [TJb)-(c)] ^2]. Note that the temperature- 
dependent Fo differs from the one determined from the 
zero voltage p^^ F = 0.210\/5 meV. A strong depen- 
dence of Pe on the gate voltage (or density) is reported 
in ^23j for parabolic two-subband wells, which also show 
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FIG. 3: Non-interacting (a)-(b) and interacting (c)-(d) Lan- 
dau fan diagrams near the v — 4 ring for several tilt angles, 
showing the anti-crossings (B & D) due to the inter-LL cou- 
pling, Eq. (|4|. For increasing Q the anti-crossings near D & 
B become larger thus making A & C move closer and the ring 
shrink. The non-interacting model [25 gives a qualitative 
picture of the ring collapse, while our many-body calculation 
(c)-(d) (T = 70 mK) provides a quantitative description and 
shows the relevance of competing XC effects, see Fig. 1(e). 



ringlike structures [23l [24] . Hence treating Fq as a fitting 
parameter is somewhat justifiable here. 

Ring shrinking and collapse. In the experiment of 
Ref. [13 the authors show the projection of the pxx side 
crossings [points A & C in Fig. [l] onto the Bj_ = B cos 
axis as a function of 0. In Fig. [TJe) we compare the 
experimental data (black dots) with the results of a non- 
interacting model (solid line) and our SDFT + Kubo 
approach (circles). The solid line is obtained from a non- 
interacting model (discussed in detail in Ref. [25^) with 
effective parameters that fit the ^ = 0° data. This simple 
model illustrates qualitatively the effects of the inter-LL 
coupling on the ring collapse, Fig. [3]^a)-(b): the anti- 
crossings [near D & B] between LLs with the same spin 
increase with 0, thus making the opposite spin LL cross- 
ings [points A & C] move closer in energy, effectively 
shrinking the ring. Note, however, that the collapsing 
angle for this non-interacting model (3.46°) is about half 
of the experimental value 0^^^ ~ 6° [13 . 

Our many-body calculation, on the other hand, agrees 
well with the data [empty circles. Fig. [TJe)] |26J. Here, 
as the tilt angle increases, a competition sets in between 
the inter-LL coupling and the exchange-enhanced spin 
splittings of the LLs. While the inter-LL coupling tends 
to shrink the rings [anticrossings. Figs. 3(a)-(b)], the 
enhanced spin splittings make the rings larger [e.g.. Figs. 
2(b)-(c) with no inter-LL coupling (i.e., = 0°)]. This 
interplay "delays" the ring collapse. Figs. 3(c)-(d); here 
we find Q^^^^'^v 6° - in agreement with [13 . 

Some remarks. Experimentally [11 all crossings of the 
V = A ring are broken at low temperature [points A, B, 
C & D] . In our approach the quantum Hall phase transi- 
tions occur only at opposite-spin LL crossings [points A 
& C]. Interestingly, we note that the inter-LL coupling 
makes the ring break near the same-spin LL crossings 
even for very small angles [see D & B in Figs. 3(c)-(d)]. 
In Ref. ^25^ we have also obtained the n2D — B map of 
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pxx for a non-interacting model and find ring breakups 
near D & B. However, the actual ring breakups near D 
& B at ^ = 0° could also be related to the derivative dis- 
continuity of the XC functionals |27j, which is absent in 
local (LSDA) and semi- local (e.g., generalized gradient 
approximations) functionals [28l [29]. This discontinu- 
ity appears as a jump in n2D at the threshold for the 
second subband occupation at 5 = [30 . The results 
of Ref. [30] suggest that orbital functionals (e.g., exact- 
exchange) [27] may give rise to phase transitions at the 
crossings B & D. Clearly more work is needed here. A 
Hartree-Fock analysis in model bilayer systems [4] sug- 
gests that quantum Hall ferromagnetic instabilities can 
occur at same-spin pseudospin LL crossings [4 . 

Summary. We have combined SDFT and linear re- 
sponse to investigate magnetotransport in 2DEGs formed 
in two-subband wells [11-13 . Our calculated n^D — B 
maps of pxx show ringlike structures. At low tempera- 
tures the V = A ring breaks due to quantum Hall ferro- 
magnetic phase transitions. The n2D — B diagram of the 
2DEG spin polarization ^ shows the ring to be ~ 50% 
spin polarized. For tilted B fields, the v = A ring shrinks 
and fully collapses at a critical angle 6>c ^ 6° , in excellent 
agreement with the data [13 . The interplay between the 
equal-spin LL anti-crossings and the XC effects are cru- 
cial here. A direct experimental evidence of our predic- 
tion of a high ^ in the ring is still lacking; the resistively- 
detected NMR data of Ref. [12 only shows signals near 
point C. We hope our work stimulates further investiga- 
tions in the literature. 
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